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ABSTRACT 

It is shown that induced Raman scattering of electromagnetic waves in the strongly 
magnetized electron-positron plasma of pulsar magnetosphere may be important for 
wave propagation and as an effective saturation mechanism for electromagnetic insta- 
bilities. The frequencies, at which strong Raman scattering occurs in the outer parts of 
magnetosphere, fall into the observed radio band. The typical threshold intensities for 
the strong Raman scattering are of the order of the observed intensities, implying that 
pulsar magnetosphere may be optically thick to Raman scattering of electromagnetic 
waves. 
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1. Introduction 

We aim to estimate the effects of the the induced Raman scattering of the strong electromagnetic 
wave propagating in pulsar magnetosphere. Raman scattering may be considered as a parametric decay 
of the initial transverse electromagnetic wave in the another electromagnetic wave and plasma wave. The 
probability of this process is greatly enhanced if there many waves present in the final states. This is 
induced Raman scattering. Another type of induced scattering - induced Brullion scattering, i.e. the 
decay of the initial transverse electromagnetic wave in the another electromagnetic wave and ion sound 
is prohibited in electron-positron plasma, since pair plasma does not support low frequency, density 
perturbing waves (like ion sound wave in electron-ion plasma) . 

Strong nonlinear coupling occurs when two waves beat together and the sum or difference frequency 
and wavelength match the frequency and wavelength of the third wave: 

W 3 = U\ + Ll>2 

k z ,3 = k z ,i + k Z:2 (1) 

In quantum language these relations may be interpreted as conservation of energy and momentum along 
the magnetic field, respectively, (the transverse components of the momentum are not conserved). The 
probability of the induced scattering depends on the number of waves in the final state, which, in turn, 
depends on the damping or escape rate of the daughter waves. The energy transfer between the modes 
will be efficient if the energy of the pump wave is strong enough to overcome the damping losses or escape 
of the generated waves. Thus, the induced Raman scattering is a threshold process: if the intensity of the 
pump wave exceeds the threshold value, the initial electromagnetic wave would start converting energy 
into the decay waves decreasing its amplitude exponentially. 
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Induced Raman scattering may be important in the pulsar environment in two ways. First, it may 
provide en effective damping of the existing electromagnetic wave, that has been generated by some 
emission mechanism at the lower altitude in the pulsar magnetosphere, where the resonant conditions 
for the induced Raman scattering were not satisfied. This may result in a short time variability which is 
generally observed in pulsar radio emission. Secondly, it may provide an effective saturation mechanism 
for the growth of the electromagnetic wave provided that the conditions for the wave excitation by some 
mechanism are satisfied in the region where an effective Raman scattering takes place. 

The first possibility, i.e. the scattering of the existent wave, is simpler to consider. As a first ap- 
proximation, we can treat the intensity of the pump wave as a constant. Then the nonlinear equations 
describing the wave coupling become linear in amplitudes of the decayed waves. The exponentially grow- 
ing solutions will imply an effective energy transfer from the pump wave. As the intensities of the decay 
waves grow, this approximation breaks down in two cases: when the amplitudes of the decay waves be- 
come comparable to the pump wave or when the amplitude of the decay waves enter nonlinear stage and 
the waves start loosing energy due to some nonlinear process (like particle trapping and acceleration). 
The net effect of any energy loss by the decay wave is the depletion of the original pump wave. For the 
practical estimates, we can neglect the nonlinear stages of evolution (like cyclic energy transfer between 
the pump and the daughter waves) and assume that if the intensity of the pump wave begins to decrease 
exponentially in the linear stage, then the wave decays completely. 

The second possibility, e.g. when the induced Raman scattering provides a nonlinear saturation mech- 
anism, is more complicated, since the intensities of all waves may be of the same order. A considerable 
simplification in this case may be obtained if the damping of one of the decay modes is very strong or 
if it leaves the region of the resonant interaction fast enough. Then, after a short period of time (when 
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the intensities of the other weakly damped waves grow considerably) the intensity of this damped mode 
is much smaller than the intensities of the two other modes and can be neglected. 

In what follows we neglect the possible nonlinear stages of the development of Langmuir turbulence 
(particle trapping or quasilinear diffusion). We also assume that the pump wave is broadband. This is 
quite different from the conventional laboratory case of a monochromatic laser-plasma interaction. The 
condition of a broadband pump wave implies that its band width Aw is much larger than the typical 
growth rate of the decay instability T: 

Acu > T (2) 

If this condition is satisfied, then we can use a random phase approximation for the statistical description 
of the interacting waves. On the other hand, the condition of weak turbulence allows one to calculate 
the matrix elements of the interaction in the approximation of stationary phases. 

The general expression for the third order nonlinear current in plasma in a magnetic field has been 



written down by [Tsytovich fc Shvartsburg 196q . An extremely complicated form of the corresponding 
expressions makes the general case of Raman scattering very difficult to consider. Several important 
simplifications can be done when considering induced Raman scattering in the pulsar magnetosphere. 
First, the superstrong magnetic field allows an expansion of the currents in 1/ujb- Second, in the pair 
plasma with the same distributions of electrons and positrons some of the nonlinear currents cancel out 
since they are proportional to the third power of the electric charge (this cancellation is exact in the 
unmagnetized electron-positron plasma). The third, less justified approximation, is that we will make is 
that all the three interacting waves propagate along magnetic field. This is an important assumption. It 
allows us to simplify the consideration considerably and to obtain some analytical estimates of the effects 
of induced Raman scattering. 
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A short overview of the work on the nonlinear process in the pulsar magnetosphere will be appropriate 
here. The possible decay processes for the transverse waves and the corresponding references are (t 
denotes transverse wave, I denotes longitudinal wave and e denotes a charged particle): (i) a decay of 
a transverse wave into two transverse wave t ->■ t' + t" flGedalin fc Machabeli 1982| ), (ii) a decay of 
a transverse wave into two Lagmiur waves t -► I + I' flGedalin fc Machabeli 19831) , (hi) a decay of a 
transverse wave into another transverse and Lagmuir wave t — > t' + I ( pcdalin fe Machabeli 1983|) Q (iv) 
induced scattering of transverse waves t + e — > t' + e QBlandford &: Scharlcmann 1976 , Sinccll <k Krolik 



19(Q, pchelkov fc Usov 1983 , Wilson 1982| ). The possible processes for Lagmuir waves are (i) a decay 



of a Lagmuir wave into another Lagmuir and transverse wave I —> t + l' QGedalin fc Machabeli 1983| ), 
(ii) a decay of a Lagmuir wave into two transverse waves I t + 1' (| Mikhailovski 198Q ), (hi) La gmuir 



wave merger / + /'—> t ( Mamradzc, Machabeli & Mclikidze 1980), (hi) induced scattering of Lagmuir 



waves ( Lyubarskii 1993| ). When treating the nonlinear process involving transverse waves, the matrix 



elements in above works have been calculated in the drift approximation - the expansion in parameter 
l/u>B was done in the very beginning (ujb = \q\B/mc is the positive nonrelativistic cyclotron frequency, 
B is magnetic field, m is mass of an electron and c is the speed of light). 

2. Kinematics of Raman Scattering in Pair Plasma 
2.1. Wave dispersion 

Normal modes of relativistic pair plasma for the case of parallel propagation consist of two transverse 
waves with the dispersion relation 

ui 2 T 

u = kc(l-5), 5 = -^ (3) 



1 The final answer for the matrix element contains probably an insignificant typographical error 
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and a plasma Langmuir wave with a dispersion relation 

- 2 = ^ + k^ft (4) 
where in Eq. (Q) the wave vectors should not be much larger than cross-over point: 

vol = k 2 Q c 2 = 2oj 2 < 7 (1 + f3f > (5) 



( Lyutikov 1998| ). We have introduced here a plasma frequency in the plasma frame uj 2 = 2a;BfJA/7 p , 



where f2 is the pulsar rotation angular velocity, A is the multiplicity factor, j p is the Lorentz factor of the 
plasma frame, T p ~< 7 > is Lorentz invariant temperature of plasma, 7 is Lorentz factor of a particle and 



fir = yTp — 1/Tp is a characteristic (dimensionless) thermal velocity of particles, is (dimensionless) 
velocity. The brackets denote averaging over the distribution function. 

At the point {uj q , k Q } the dispersion curve of plasma waves intersect that of waves in vacuum, so 
that for k > k D the plasma waves become subluminous. The necessary requirement for that is that the 
distribution function falls fast enough at large values of momentum (faster than 1/p 4 ) ( [Lominadze & 
Mi khailovskii 1978| ). 



In the subluminous region, the plasma wave dispersion relation may be written as 

uj = kc — rj{k — k )c, where 



77 = 1 - 

2.2. Scattered frequencies 

From the resonant conditions 



de zz \ lfde z 



dk J / V da; 



UJo,ko 



<7(i + /?r> _L 
<7 3 (i+/?) 3 > ~ n 



(6) 



k 3 = h + k 2 (7) 
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it follows that that the only kinematically allowed Raman scattering process is a back scattering (k\ < 0) 
of the initial electromagnetic wave (see Fig. 1). We find then 

k3 = Yc{r^ + k2C 



2c V 1 - 6 

In the region just above the cross-over point {uj , k Q } it is possible to use the expansion (|6|) for the 
dispersion of the plasma waves. The resonant wave vectors are then given by 

, k 2 + k / k 2 -k a \ k 2 -k Q 

The wave vectors in Eq. (^) are limited by k 2 — k Q < k D , \k\\ < k Q . 

2.3. Induced Raman Versus Induced Compton Scattering 

Induced Raman scattering is closely related to the induced Compton scattering of waves by the plasma 



particles. In the conventional treatment of induced Compton scattering in pulsar setting (e.g. [Blandford 



fe Scharlemann 1976 , Sincell fc Krolik 1992 ) the collective effects of the plasma are ignored. In the case of 
induced Raman scattering collective effects play a major role. Raman scattering is a limit of the induced 
Compton scattering of waves by particles in a medium when the beat wave of the incoming and scattered 
wave becomes a normal mode of the medium (equivalently, it falls on the "mass surface" of the medium 
dispersion equation). 

In the nonrelativistic plasma the collective treatment of the wave scattering by plasma particles is 
justified if the following condition is true: 

Ak\ D < 1 (10) 

where A A; is the change in the wave number of the scattered wave and Ad is the Debye length. Condition 
( p!o[ ) states that the wave number of the oscillations of the electron on the beat of the two scattered 



waves be much less than the inverse of the Debye length. If condition (|l^) is not satisfied, then the beat 
wave will either experience a strong Cherenkov damping ( if AkXp ~ 1) or the beat wave will not feel 
the presence of a medium ( if AkXo 3> 1), so that the scattering process will be described as induced 
Compton scattering. 

In the relativistic plasma with the characteristic Lorentz factor of the particles < 7 > 3> 1 the condition 
( |l0| ) is modified. To find the the generalization of condition (1C) to the case of relativistic plasma we 



note, that condition ( JlOj ) may be rewritten in therm of the phase velocities of the beat wave v p h = uj/k 
(which is a plasma wave in case of Raman scattering) and the thermal velocities of the plasma particles 
vt =< v 2 > 1 / 2 . Condition (|lC|) then can be rewritten 

Vph » v T (11) 

Using (||), we find that the phase velocity of the subluminous plasma waves becomes equal to the 
thermal velocity of plasma particles approximately at k ~ 2k Q . Thus, the following condition should hold 
for the wave-particles scattering in relativistic pair plasma to be considered as a Raman scattering: 



Ak < 2k Q ps V2 < 7 >uj p (12) 

When Ak ~ 2k Q the plasma waves are strongly damped and when Ak ^> 2k Q the wave-particles scattering 
is described as single particle Compton scattering. As we will see later (Eq. (|4l|)) the condition ( |l2|) is 
satisfied in our problem. 

3. Transition probabilities 

First we calculate the matrix elements for the nonlinear three wave interaction in the strongly mag- 
netized electron-positron plasma when all three waves have their wave vectors along the magnetic field. 



We represent the electric fields of the the three interacting waves in the form 

E(r, t) = U [Ei exp{i(k ir - u x t)} + E 2 exp{i(k 2 r - u 2 t)} + E 3 exp{i(k 3 r - uj 3 t)}} (13) 

The indices 1,2,3 refer to the scattered electromagnetic wave, plasma wave and the initial electromag- 
netic wave respectively, K denotes the real part. We consider lineally polarized waves. For the parallel 
propagation the case of circularly polarized waves may be trivially reduced to the scattering of two lin- 
eally polarized waves. Implicit assumption here is that the coupling between the waves is weak, so that 
the waves retain their identity as eigenmodes of the medium during the times much longer than inverse 
of their frequencies. 

We calculate transition probabilities by finding the coupling coefficients for the resonant three-wave 
interaction. The interaction of two transverse and one longitudinal wave propagating along magnetic field 



in a relativistic plasma has been considered by Brodin &z Stenflo 1989 . For the decay process uj 3 = u\ + u> 2 
and k 3 = k\ + k 2 the coupling of the waves is described by the system of equations 



dt 
dE 2 

dt 
dE 3 

dt 

where 



c x E* 2 E 3 
c 2 E\ E 3 

c 3 E{ E 3 (14) 



^1 r 



C2 = -yr — 



and 
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dF(p) 



V J 
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(oj - kv z )p± dF(p) 

uj — kv z =F u) c dp± ' 7717(0; — kv z =F u> c ) dp z 
m dF(p) 
7^2 dp z 



=F 







dp 2 . 



UcP± 
7 4 W2 



fcl 



27 3 mc 2 cD2 

^3 



+ 



(&1 =F w c ) (u) 3 =F w c ) (cDi =F w c ) (w 3 =F w c ) 
7 2 (wi =F uj c ) (lu 3 =f w c ) 



uj r h 



c^3 



^3 =F w c ) (cji T w c ) (u>i =F u; c ) a> 2 (& 3 =F u c ) uj 2 



C 2 



ap(p)i 
J 



*Xp) 



(16) 



where 0)1,2,3 = ^1,2,3 — ^1,2, 3^z an d w c = qB /^mc is the relativistic cyclotron frequency. 

We now simplify this relations for the case of one dimensional pair plasma with the same distribution 
functions. Setting F(p) = 2 \S(p 2 L )f(p z ) and introducing a positive nonrelativistic gyrofrequency u>b = 
a q qBo/mc (a q is the sign of the charge) we find 



1 + 



H f dp z 8f(p z ) 



k 



1 - 2uj1 



u> — kv z dp z 
dp z (uj -kv z ) 2 
7 (u> - kv z ) 2 - L0 2 B /~f 2 



f(Pz 



c 



2uj p \q\ f dp z lo b 



uoiuo 3 m 



k 3 Cj x 



k\(b 3 



7 5 (b\ \ cJi 2 



f(Pz 



(17) 
(18) 
(19) 



The probability of the decay is ( [Melrose 1978c] Eq. 10.124) 



7^1 E>2 E>3 

n(k 3 , k 2 , ki) = (27r) 5 ft E E E Icui^sCI 2 5(k 3 - k 2 - ki) 5(co 3 -u 2 - u x ) 



ijJ 3 UJ\UJ 2 

dp z ( k 3 toi 



(2Tr) 5 ?iq 2 u;y B u; 2 
2m 2 u>iuj 3 

5(k 3 - k 2 - ki) 5(uj 3 -u 2 — U\ 



kiu 3 



7 5 u>| ^cJi 2 - wl/7 2 



^3 2 - w|/ 7 2 



fiPz 



(20) 
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with R\ = R% = R% = 1/2 
In the cold case we have 



n(k 3 ,k 2 ,ki) 



(2vr) 5 V^H 
2m 2 (jJiUJ2,uJ2 



5(k 3 - k 2 - ki) 5(uj 3 -L02- uj\) 



(21) 



which may be obtained from Melrose 1978c Eq. 10.105. 

Equations ^ and (H) can be simplified in the limit of a very strong magnetic field. In the limit 



2ui 



OJ 



B 



1 H ^ / dp z ^(uj - kv z ) 2 f{p z ) 



C = ±- 



2u;l\q\ f dp z 1 fk 3 h 



rn 



7 3 Qj\ijJb V^3 ^1 



f(Pz 



(22) 
(23) 



Using Eq. ([l?]) the equation (^3| ) may be rewritten in the form 



For the back scattering Eq. (23) gives 



C 



, 2|g|fc 2 
± (e 2 



(24) 



(25) 



For the plasma waves on the mass surface e zz = and 



C = ±- 



2|^|fc 2 



n(k 3 ,k 2 ,ki 



2 2 g 2 A:|(27r) 5 ^ 



2,. ,2 



u; 2 <5(k3 — k 2 — ki) (5(u;3 - o» 2 — cji) 



(26) 



We can compare this probability with the unmagnetized case ( Tsytovich 1970| ). In the strongly mag- 
netized plasma the probability of Raman scattering is decreased by a ratio uj 2 3 /uj b . This is similar to 
the suppression of the Thompson scattering, whose probability is decreased by the same ratio. 



11 



It is convenient to describe the wave distribution in terms of photon occupation numbers 

F 2 

" k = Rt) (27 > 

where is the spectral energy density, which is related to the total energy density 

w = < 28 > 



The kinetic equation for the Raman scattering are (Melrose 1978c, Eq. (5.93)) 

^^ = " /tSf/ 7 |^3"(k3,k 1 ,k 2 )K(k 1 )n 2 (k 2 )-n3(k) (n 1 (k 1 )+n 2 (k 2 ))) 

+r 3 n 3 (k 3 ), (29) 

^aT^ = ~ IW^I (^^ k 3' k i' k 2)(^( k i)^(k 2 )-n 3 (k3)(n 1 (k 1 )+n 2 (k 2 ))) 

+r 2 n 2 (k 2 ), (30) 
f^sf 7 ^y3n(k3,k 1 ,k 2 )(n 1 (k 1 )n 2 (k 2 )-n 3 (k3)(n 1 (k 1 )+n 2 (k 2 ))) 

+r 3 n 3 (k 3 ) (31) 

Where I\ 2) 3 are the corresponding damping or growth rates. The most unimportant damping is the 
Landau damping of the plasma waves n 2 (k 2 ) on the particles of plasma. In ( |3l| ) the total derivative on 
the left hand side is 

d d duj d duj d 

dt = di + dkd^~'d^dk ^ ' 

with the second term describing the convective transport of the waves. 

Next we integrate over transverse wave vectors to find the system of equations describing the scattering 
of the waves propagating along magnetic field 



dn\{fa 
d~T 



dfa f dfa 



(2tt)J (2tt) 



u(fa, fa, fa) (ni(/ci)n 2 (/c 2 ) - n 3 (fa) (n^fa) + n 2 (fa))) 
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where 



and 



+rmi(fci), (33) 

dn2 ^ = ~ J J J^j u ( k2 > k ^ ( n i( fc i) n 2(A;2) - n 3 (k 3 ) (ni(fci) + n 2 (fc 2 ))) 

+r 2 n 2 (/c 2 ), (34) 
dn 3 (k 3 ) /" dfei /" dk 2 



(jf J ,9-)/ ( -, )7r y u ( fc 2, fei, fe) («i(A;i)n 2 (fc 2 ) - n 3 (& 3 ) (ni(fei) + n 2 {k 2 ))) 

+T 3 n 3 (k 3 ) (35) 



n 



(h) = J ^^ ni (k t ), i = 1,2,3 (36) 



2Vfc 2 2 (2^) 5 ^ 2 ^ 

u(fc 2 , fci, fc 2 ) = 5—3 o(fe 3 - fc 2 - fei) d(u;3 - a> 2 - cji) (37) 

The two (5- functions in u(k 2 , fei, A; 2 ) allow us to perform integrations in (]3 



dmteO _ 2^^1(2^)3^ (ni(Afi)na(Jfca) _ n3{h) {ni{ki) + Mh))) + r.mih), (38) 



dt m 2 uj\ 
dn 2 (k 2 ) 2 2 q 2 kl(2-Kfhuj 2 



fdun\ 



(ni(ki)n 2 (k 2 ) - n 3 (k 3 ) (ni(fci) + n 2 (k 2 ))) + T 2 n 2 (k 2 ), (39) 



dn 3 (/c 3 ) _ 2 2 q 2 kj(2ir) 3 huj 



dt m 2 io 2 



B 



( §^2_\ 

\dk 2 ) 



2 {ni(k 1 )n 2 (k 2 ) - n 3 (k 3 ) (m(ki) + n 2 {k 2 ))) + T 3 n 3 (k 3 ) (40) 



Implicit assumption in Eqs. ( |3Sj ) - (|40| ) is that the wave vectors on the right hand side of Eqs. (|38|), 
(|39|) and (|io|) satisfy the resonance conditions. Thus in Eq. (|38|), for example, the two independent 
variables are time t and wave vector k±. The explicit relations between wave vectors, corresponding to 
Eqs. (H), d|) and ©, are 

k 3 ^k 2 = 2T 2 \k 1 \ + k Q , uj 2 = (2T 2 \k 1 \ +k - 2|fci|)c 
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,,, k 2 k . k 2 k Q 
\ki\ = 2J , 2 , fc3~K 2 , uJ 2 = {k 2 ^— )c 

I fell ~ k3 2T 2° > k 2 ^k 3 , u 2 = (k 3 - k3 T2 k ° )c (41) 



respectively for Eqs. (pq), (|39| ) and (40). The quantity /c Q = y2Tu) p is the cross-over wave vector, where 
the dispersion relation for the plasma waves intersect the vacuum dispersion relation. 

From Eqs. and (JT2|) we find that the induced Raman scattering occurs in the frequency region 



Vzfujp <kc< 2sf 7 2Tu) v (42) 
This frequency window (in the plasma frame) corresponds to the frequencies in the pulsar frame in the 



range 2^JTXujb^I/j p < kc < A^jTXojB^/lp-, which is 5 x 10 8 radsec" - 1 < kc < 10 9 rad sec which falls 
into the observed frequency range. 

4. Application to Pulsars 

In this section we consider the possible applications of the induced Raman scattering to the wave 
propagation and spectrum formation in the pulsar magnetospheres. 

4.1. Scattering of the Pump Wave 

There is threshold intensity of the pump wave that defines what is known as strong induced Raman 
scattering. The rate of induces plasmon production, which depends on the intensity of the pump wave, 
should exceed the rate of loss due to the various damping processes. The end result of such strong 
scattering is a total depletion of the initial pump wave. The threshold photon density for strong Raman 
scattering is determined from ( |39[) 

dki f dks 



■u(k3,ki,k 2 )n3(k 3 ) =T 2 (43) 

Z7T J Z7T 
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where r 2 is the damping rate of plasma waves. 

Linear Landau damping rate in the relativistic plasma (T p 3> 1) is ( [Lyutikov 199S 



7r<4 / 3 df 



Tkcuj 2 V #7 

with 



(44) 



h Q T 2 

7 ^ es = ~AlP k °° = ^ = V 2TpWp ' Ak = k-k (45) 

The equations are valid for A/c < fco- 

With a power law distribution function / oc j~ a , a > 2 we find 

(I t>2\ -a/2+1 
^) W 

Note that as Ak — ► 0, i.e. when the wave vector approaches the cross-over point where the phase speed 
is equal to the speed of light, r 2 ^ since there no particles with the speed equal to the speed of light. 
The threshold spectrum of the incoming radiation is then given by 

Estimating Ak ~ k Q we find the characteristic threshold photon density: 

nf - 10- 2 ^^ T- 1 - = 2.7 x 10 22 ^ £ 12 A 5 ^ P . 2 Tf 3 i? 9 3 7p , 10 (48) 

where S i2 = B/W 12 G is a magnetic field at the surface of the neutron star, A5 = A/10 5 is the multiplicity 
factor, i.e. the ratio of the plasma density to the Goldreich- Julian density uqj = fi • B/(2 7rec), -P0.2 = 
P/0.2sec is the pulsar period, T\q = T/W is the relativistic invariant temperature, Rg = i?/10 9 cm is the 
distance from the neutron star, 7p,io = 7p/10 is the average streaming Lorentz factor of the secondary 
plasma and we used a = 2. Note, that this is an upper bound on the threshold photon density (for 
k ~ 2k Q ). At the lower frequencies k ~ k Q the threshold values are much lower. 
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For the relativistic Gaussian-type distribution / oc 1/T x exp{— 7/T} the damping rate is 



(compare with Lominadze fc Mikhailovskii 1978 Eq. (2.10)). The spectrum of the incoming radiation is 
then given by 



/Afcr 3/2 [ /2^T^ 1/4 ' 
128vrV^ 9 / 4 ^ ^ Up J exp | ^Ak 2 , 



n 3 « ,ooll!L ( — ) exp <( - ( ) }> (50) 



The characteristic threshold photon density is in this case 

232 23 

„th ^ ln -4 m C U B -in~5 m C u Blp _ 9 7 v ln 20 1 R x -l p T -3 p-3„, / M \ 

3 ~ q 2 T 3 hu 2 ~ q 2 T 3 hQX ~ 12 5 a2 10 9 7p ' 10 ( ^ 

(upper bound) where we used u> 2 = 2lubQX/j p - 

For the comparable values of the temperature, i.e. for the similar dispersion of the energies of the 
plasma particles, the threshold in the case of a power law distribution is much larger than in the case of 
Maxwell-type distribution. This is due to the larger damping rate in the case of a power law distribution 
since there are more particles satisfying the Cherenkov resonance in this case than in the case of Maxwell- 
type distribution. 



The estimates (48) and (51) are done for the one dimensional (in k) photon density in the plasma 

frame. The photon densities ( |48| ) and (|5l|) are the one dimensional photon densities. 

The energy flux per unit interval of frequencies is given in terms of photon density as follows: 

n(k,r)chujdk n(k,r)chujdk 

b [y)dv = = = = n{k, r)nojdv (52) 

Sir 6 2tt 

Where we took into account that the dispersion relation of the electromagnetic waves is almost vacuum 
like, so k ~ co/c. 
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We recall next, that both three dimensional occupation numbers n(k, r) and one dimensional photon 
density n(k,r) (Eq. (refkd921)) are relativistic invariants. It follows from (|52| ) that 



F{v) 



mv 



(53) 



(for the one dimensional flux, i.e. integrated over transverse wave vectors). Thus, the flux in the observer 
frame is 



F'{y') = 2im(k',r)hv' in 



erg 



cm 2 sec Hz 



We can estimate the observed flux (in Jankys) at the Earth: 



F obs (i/) = F'(u" (R 



Using Eqn. ©, (||), (||) and (||) we find 



F 



obs 



1Q- 2 



1CT 4 



n;m z c 6 u>* B 'ypV f Rns \ ( Rns \ 
q 2 n\ V d J V R J 



for the power law and exponential distributions correspondingly. 



nobs 



V 



Jy 



(54) 



(55) 



(56) 



(57) 



Numerical estimates give 

/ 

30 

3 x kt 1 

These are comparable with the fluxes observed from bright pulsars. Since the estimated threshold in- 
tensities are upper bounds, we conclude that stimulated Raman scattering may be an important factor 
in the formation of pulsar radio emission. Estimates for the millisecond pulsars give about an order of 
magnitude lower values. 
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4.2. Saturation of Cyclotron Instability by the Induced Raman Scattering 

Cyclotron-Cherenkov instability, operating in the outer parts of magnetosphere at the anomalous 
Doppler resonance 

u(k) - k\\v\\ - s— = for s < (58) 

7 



has been suggested as a possible mechanism of pulsar radio emission ( [Kawamura &: Suzuki 1977 , Lomi- 



nadze, Machabeli & Mikhailovskii 1979, Kazbegi et al. 1991b, Lyutikov 1998). In this section we consider 



the nonlinear saturation of the cyclotron-Cherenkov instability by the induced Raman Scattering. To 
simplify the consideration, we assumes that the product electromagnetic waves (labeled by index 1) leave 
the region of resonant interaction, so we can assume that their photon density is ni(k\) ~ 0. We are 
interested in the time asymptotic limit for the intensities of the transverse waves, excited by some insta- 



bility mechanism with a growth rate IV Neglecting the wave convection we then find from Eqs. (38) 
and ©: 



2 2 q 2 B(2ir) 3 hu 2 „ N „ „ , 
m ^(L A 2 ns(k 3)+ T 2 n 2 (k 2 ) = 0, (59) 

2 2 q 2 k 2 l (27r) 3 huj 2 . . 
171 ^ B {dk^J 

Comparing Eqs. ( f43| ) and (^) we find that they are identical: the threshold intensity for the strong 
Raman scattering is equal to the saturation intensity. This is not a surprising fact, since in both cases 



it is the damping of the plasma waves that controls the electromagnetic wave intensity. In Section [4J. 
we found that the corresponding intensities are comparable to the observed ones, so that the nonlinear 
saturation by the induced Raman scattering may provide an observed emissivities. 

We also note in this context, that at a given radius, the induced Raman scattering operates in a lim- 
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ited frequency range (12). So, if the instability that produces radio emission is intrinsically broad band 
(with the bandwidth of the growing mode larger than the range in which the induced Raman scattering 
operates), then the effects of the induced Raman scattering would produce a dip in the spectrum, corre- 
sponding to that frequency range. If, on the other hand, the instability that produces radio emission is 
intrinsically narrow band with the central frequency of the instability falling into the operational range of 
the induced Raman scattering and changing with radius to mimic the observed broadband spectrum, then 
the induced Raman scattering can be an effective saturation mechanism for a larger range of frequencies. 

5. Conclusion 

In this work we considered induced Raman scattering of a transverse electromagnetic wave propagating 
along magnetic field in a pair plasma. We found that induced Raman scattering can be an important 
factor for the wave propagation in the pulsar magnetosphere and as a nonlinear saturation mechanism for 
the electromagnetic instabilities ( cyclotron- Cherenkov instability, for example). The threshold intensities 
for strong Raman scattering, which are determined by the damping rate of the Langmuir waves in the 
relativistic plasma, are comparable to the observed fluxes from radio pulsars. 

One of the major limitations of this work is the assumption of the propagation along the magnetic field. 
The natural extension of this work to oblique propagation would involve several complications: (i) there 
will be a distinction in the transition probabilities between extraordinary mode (with the electric field of 
the waves perpendicular to the k — B plane) and ordinary mode (with the electric field in the k— B plane), 
(ii) ordinary mode has an sensitive dependence on the angle of propagation, so that for the ordinary mode 
the Raman scattering will occur only in a very limited angle around magnetic field; the scattered waves 
will be a backward propagating, (iii) for the extraordinary mode the Raman scattering is possible for 
large angles of propagation; in this case both electromagnetic waves can be forward propagating. 
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FIGURE CAPTIONS. 

Fig. 1. Kinematics of backwards Raman scattering in the plasma frame. The vacuum dispersion 
relation uo = kc, the transverse wave dispersion relations u> = kc(l — 5) and the plasma wave dispersion 

/ q \ 1/2 

relations (which starts from ( ^ J u) p ) for the parallel propagation are shown. The bold faced vectors 
represent k = {uj, k}. Cross-over point is given by k Q and u . Only waves with k > k t can be scattered. 
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